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Abstract 

A phenomenological characterization of hysteresis in magnetostrictive materials is pre- 
sented. Such hysteresis is due to both the driving magnetic fields and stress relations within 
the material and is significant throughout most of the drive range of magnetostrictive trans- 
ducers. An accurate characterization of the hysteresis and material nonlinearities is necessary 
to fully utilize the actuator /sensor capabilities of the magnetostrictive materials. Such a 
characterization is made here in the context of generalized Preisach operators. This yields 
a framework amenable to proving the well-posedness of structural models that incorporate 
the magnetostrictive transducers. It also provides a natural setting in which to develop prac- 
tical approximation techniques. An example illustrating this framework in the context of a 
Timoshenko beam model is presented. 
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1 Introduction 


Important members in the class of smart materials currently employed in control applica- 
tions are magnetostrictive materials. The phenomenon of magnetostriction is defined as the 
strain which results when a magnetic material is subjected to a magnetic field. While this 
phenomenon occurs in most ferromagnetic materials, only in recently developed rare-earth 
materials are the strains and forces sufficiently large to facilitate their use in actuators and 
sensors. This has led to the use of magnetostrictive materials as ultrasonic transducers, posi- 
tioners, sonar projectors (500-2000 Hz) and isolators (5-60 Hz). They are also being considered 
as actuators for controlling vibrations in thick structures and in heavy rotating components 
such as milling machine bits (further details regarding applications in which magnetostrictive 
transducers are employed as sensors and actuators can be found in [8, 16]). 

To fully utilize their capabilities as either sensors or actuators, the input and output 
characteristics of magnetostrictive materials must be quantified in a manner amenable to 
parameter estimation and control applications. As detailed in the literature (e.g., see [17, 19]) 
and outlined in the next section, the input/output characteristics of magnetostrictive materials 
are inherently nonlinear arid display significant hysteresis. To utilize the transducer responses 
in a range useful for full control applications, accurate mathematical characterization of the 
nonlinearities and hysteresis is required. 

The general techniques we use here to characterize the magnetic nonlinearities and hys- 
teresis are based on Preisach-type operators. The use of Preisach techniques for characterizing 
general magnetic hysteresis is well-established in the literature (see [1, 6, 13, 15, 23, 24, 26] 
and references therein) and some aspects have been extended to magnetostrictives [2, 9]. In 
[2], the classical Preisach model w r a.s modified for Terfenol-D materials by considering a char- 
acterization for the output strain in terms of two inputs and two Preisach kernels. This was 
motivated by the observation that magnetostrictive materials exhibit hysteresis with respect 
to both magnetic field and stress. Further coupling between the kernels and input was included 
in the model proposed in [9]. 

In this work, we consider the modeling of nonlinearities and hysteresis in magnetostrictive 
transducers in the context of the generalized Preisach or Krasnoselskii-Pokrovskii kernels 
developed in [6]. Such generalizations of the classical Preisach theory are motivated by the 
goal of attaining a kernel which is continuous with respect to both time and shape parameters 
(as proven in [6], the classical Preisach kernel is discontinuous in both aspects). This provides 
a framework in which to prove well-posedness of models which incorporate magnetostrictive 
transducers. Furthermore, it yields a framework which is natural for the development of 
approximation techniques for both simulations and parameter estimation. 

It should be noted that the modeling of nonlinearities and hysteresis through Preisach 
techniques is phenomenological rather than physics-based. It provides a mathematical charac- 
terization of the input/output relationships for magnetostrictive materials in lieu of a complete 
theory for the electromagnetic and magnetomechanical properties of the materials. As the 
physical theory is advanced, physics-based models may provide additional insights regarding 
the actuator/sensor mechanisms and hence replace the Preisach models. For example, theory 
based upon magnetic domain wall motion, which is used in [21, 22] to model hysteresis in fer- 
romagnetic materials, may be applicable to magnetostrictive materials. At the current time. 
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however, several physical mechanisms in magnetostrictive materials are still not well under- 
stood, thus motivating the phenomenological or empirical characterization of the material. 

The construction of magnetostrictive transducers and physical properties of the magne- 
tostrictive materials are described in Section 2. The source and form of the nonlinearities 
and hysteresis are detailed to illustrate issues to be addressed in the models. The generalized 
Preisach kernels and operators developed in [6] are summarized in Section 3. In Section 4, 
a thick cantilever beam with magnetostrictive actuators is used as a prototype for structural 
models which incorporate magnetostrictive inputs. Modified Timoshenko equations are used 
to model the beam dynamics while generalized Preisach or Krasnoselskii-Pokrovskii opera- 
tors are used to characterize the inputs. The well-posedness of the model and parameter 
estimation problem is established using the framework of [4, 6]. Although the framework is 
described here in the context of a specific, simple structure, it is sufficiently general to include 
structures comprised of plates and shells as well as structures in which nonlinear dynamics 
are significant. 

2 Magnetostrictive Transducers 

To illustrate issues which must be addressed when modeling magnetostrictive materials, the 
transducer depicted in Figure 1 is considered. This transducer is typical of those currently 
employed in applications and contains the basic components required for experimental use as 
a control actuator. Details regarding the construction and performance of this transducer can 
be found in [17]. 

The primary components of the transducer consist of a magnetostrictive rod, a wound wire 
solenoid, and a cylindrical permanent magnet. In current transducers, the magnetostrictive 
material is typically composed of terbium and dysprosium alloyed with iron. A commonly 
employed material is Terfenol-D (Ter: terbium, fe: iron, nol: Naval Ordinance Laboratory, D: 
dysprosium) which is constructed as a cylindrical rod and placed in the center of the trans- 
ducer. The sensor/actuator capabilities of the material are due to the rotation of magnetic 
moments within the rod in the presence of an applied magnetic field. In ferromagnetic ma- 
terials such as Terfenol-D, moments are ordered and exhibit a high degree of alignment at 
temperatures below the Curie point (regions in which moments are highly aligned are termed 
domains). As depicted in Figure 2a, the moments within domains are primarily oriented per- 
pendicular to the longitudinal rod in the absence of an applied field. Prestressing the rod with 
the spring washer serves to increase the number of moments perpendicular to the axis (see 
Figure 2b) and places the material in compression. This latter objective is necessary due to 
the inherent brittleness of Terfenol-D. In the presence of a magnetic filed, the moments rotate 
so as to align with the field. Consequently, if the field is applied in the direction of the rod 
axis, the moments rotate in the sense depicted in Figure 2c and significant strains are gener- 
ated. This is termed the Joule effect and provides the actuator capabilities of the transducer. 
Sensing is accomplished through the measurement of the magnetic fields which result when 
mechanical stresses cause rotations of the moments within domains (Villari effect). Details 
regarding these effects can be found in [19, 25]. 

The strains generated through an applied field are always positive since rotation of the 
moments from the prestressed perpendicular state leads to an increase in the rod length. As 
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Wound Wire Solenoid 


Figure 1 . Cross section of a typical Terfenol-D rnagnetostrictive transducer. 

indicated in Figure 3, the relationship between the applied magnetic flux or induction B and 
strain t is also highly nonlinear with saturation occurring at large field strengths. Moreover, 
slight hysteresis also exists between B and e at high drive levels (this is not depicted in the 
figure). 

The generation of bidirectional strains is accomplished through either a DC current 7 0 
applied to the solenoid which surrounds the rod, or an enclosing cylindrical permanent magnet 
which provides a biasing magnetic induction Bo. A time varying current I(t) is then used to 
vary the induction in the rod between 0 and B m . This provides the capability of generating 
both positive and negative strains. 



Figure 2. Magnetic domains in the Terfenol-D rod; (a) Orientation of moments within 
domains in unstressed rod with no applied field; (b) Orientation of moments in prestressed 
rod with no applied field; (c) Orientation of moments in prestressed rod when field is applied 
in direction of longitudinal rod axis. 
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Figure 3. Strain distribution t generated by an applied magnetic induction B. 

To model the transducer for actuator and sensor purposes, it is necessary to characterize 
the relationship between the current I applied to the solenoid, the resulting magnetic field //, 
the associated magnetic induction B and finally, the generated strains e. Furthermore, the 
quantification must incorporate the contributions due to the permanent magnet. 

As detailed in [19], the magnetic induction and field are related by the permeability // 

which is defined as 

B 



The magnetic field is due both to the solenoid and the permanent magnet. The magnetic 
induction from the permanent magnet is approximated by B 0 = while Ampere’s law 

yields B = fml, where n is the number of turns per unit length in the solenoid, as the 
approximate magnetic induction due to the solenoid (these are approximate since edge effects, 
air gaps, etc., are neglected). 

For Terfenol-D, the permeability fi is highly nonlinear and exhibits significant hysteresis 
as indicated by the induction/field relations depicted in Figure 4a. As discussed in [2], this 
hysteresis is manifested with respect to both the applied magnetic field and stresses within 
the magnetostrictive material. To indicate the latter contribution, the permeability in mag- 
netostrictive applications is often denoted by /.r . The hvsteretic relationship between the 
magnetic field, magnetic induction and material stress are then inherently manifested in the 
field-strain relations as shown in Figure 4b. 

As indicated by the preceding discussion, the nonlinear relationships between the applied 
current and magnetic induction are augmented by nonlinearities in the behavior of the magne- 
tostrictive materials. For example, experimental results in [11, 14] indicate that the Young’s 
modulus E h for Terfenol is dependent upon the applied magnetic field which partially accounts 
for the dependence of magnetic hysteresis on the material stress. Furthermore, experimen- 
tal results in [11, 14] demonstrate that other material properties such as magnetomechanical 
coupling coefficients are highly sensitive to operating conditions such as prestress level, AC 
drive levels, operating frequencies and temperature. To provide a framework amenable to 
characterization of these relationships, we summarize next, necessary theory regarding gener- 
alized Preisach operators. The application of this theory to magnetostrictive materials in a 
structural setting is considered in the final section. 
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(a) (b) 


Figure 4. (a) Relationship between the magnetic field strength H and the magnetic flux 
density B\ (b) Applied magnetic field H and resulting strain distribution e. 

3 Preisach and Krasnoselskii-Pokrovskii Operators 

In this section, we provide basic theory regarding Preisach-type operators in a Hilbert space 
setting. Classical Preisach kernels and operators are defined first and are then extended to 
kernels of Krasnoselskii-Pokrovskii type to attain desired continuity properties. The material 
in this section summarizes theory from [6, 7], and the reader is referred to those references for 
further details. 

To motivate the general kernels used later, we first illustrate with a single delayed, relay 
operator k. This kernel is characterized in terms of crossing times r(t) defined by 

T(t) = {?/ € (0, T] I u(tj) = Si or u{i]) = s 2 } 

where s = (si,.s 2 ) are points in the Preisach half plane 

<S = p G IR 2 | s = ($i, s 2 ), < s 2 | 

and u denotes an input function. The values si,.s 2 are threshold values for the multivalued 
kernel as reflected in the definition 

[M«>0]( 0) if t(<) = 0 

[fc s (u, £)](*) = < — 1 if r(f) ^ 0 and «(raaxr(f)) = si 

k +1 if r(t) 7^ 0 and u(maxr(f)) = s 2 . 

A depiction of this kernel is given in Figure 5. The starting value 

{ —1 if u(0) < Si 

if if s i < u (0) < ,s 2 

+ 1 if «(0) > »2 

defines the initial state of the kernel in terms of the parameter £ £ { — 1, 1}. 
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The output remains on a branch until a threshold is reached in the monotonically increasing 
input u. At that point, the output jumps to the other saturation value and remains there 
until the other threshold value is reached. For example, an output response starting with a 
value of —1 will retain that, value until u(t) reaches s 2 . The output then jumps to +1 until 
the threshold value of Si is reached. 

The classical Preisach operators are then defined in terms of parallel collections of these 
single relay operators. To this end, we let M. denote the set of all finite, signed Borel measures 
on S and let / be a Borel measurable function mapping <5 — > { — 1,1}. For u € C[0, T) and 
// £ j\4, the Preisach operator is defined by 

[P,(ti,/)]( t) = J s [k,(u,f{s))]{t)dv(»). 

The goal in the parameter identification problem is to estimate the measure v so that the 
model response “fits'' experimental data in a least squares sense. 

While this provides an operator which is useful for many applications, this classical defi- 
nition does not yield a kernel, and hence operator, which is continuous with respect to either 
time or parameters. Specifically, as proven in [6], the mapping in time 

t [k s (u,£)}{t) 


and the parameter space mapping 

s i y [k s (u, £)](f ) 

are discontinuous for the classical Preisach kernel k s (u,£). Continuity in time is important 
from a physical perspective while continuous parameter dependence is crucial for the devel- 
opment of practical parameter estimation methods. 



Figure 5. Single Preisach relay operator with threshold values .si,,s 2 . 
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To avoid the difficulties associated with the discontinuous mappings, a Krasnoselskii- 
Pokrovskii kernel of the type discussed in [6] is employed. This kernel is somewhat less 
general than the influence operators considered in [23] and arises as an extension of smoothed 
Preisach operators. These operators differ from the previously-defined Preisach operator in 
the manner through which an envelope of admissible paths is defined. In this case, an envelope 
is provided by translates 

r Sl = r(x - sj) 
r S2 = r(x - .s 2 ) 

of a Lipschitz continuous ridge function r(;r) as depicted in Figure 6. For monotone inputs 
u m € C f [0, T], a monotone output operator is defined by 


[7£(u m , £)](*) 


max{(, r(u m (t) — s 2 )} if u m is non-decreasing 
min{£, r(i/ m (/) — S\)} if u m is non-increasing. 


In terms of this operator, a kernel is defined for piecewise monotone inputs ti pm G C f [0, T] fl 
5'ij[0, T], where ^ij[0,T] is the set of piecewise linear splines with j knots in [0, T], in the 
following manner. The initial value of the operator is taken to be 7Zo = £. A kernel k s is then 
defined recursively on each subinterval by 


[k,{Up m X)](t) = [R-(Upm,'R'k~l)]{t) , t e [4-1,4] (3.1) 

where 7 Zk = TZ(u pm , lZk - 1 )(ik ) , k = 1, * * * ,j. The input and action of this kernel are illustrated 
in Figure 7. This provides a definition of the kernel useful for computational algorithms in 
which inputs are discretized in terms of a piecewise linear basis. This definition is readily 
extended to arbitrary a G F f [0,T] through standard density arguments as detailed in [23]. 

It is natural in the applications considered here to formulate the hvsteretic input operator 
in a manner commensurate with a weak model formulation. We thus consider a space of test 
functions V and state space H in which V is continuously and densely embedded in H and 
forms a Gelfand triple 

V H ~ H* ^ V*. 



Figure 6* Hysteresis envelope provided bv the translates r\ and V 2 of the ridge function r(x). 



(a) (b) 


Figure 7. (a) Piecewise monotone input; (b) Output from the Krasnoselskii-Pokrovskii kernel 
in response to a, piecewise monotone input. 

The following theorem from [6] summarizes the continuity properties of the kernel k s and 
quantifies the resulting input operator B u . 

Theorem 1. Let k s denote the kernel defined in (3.1) and let S denote the closure of S in IR 2 . 

(1) For each s e S and £ G {—l, 1}, the kernel satisfies k s (-,t ;) : C’[0,T] — ► C[0, T] . 

(2) For each u G C[0,T],£ € {-1,1} and t € C[Q,T], the map s [ k s {u,£)){t ) is continuous 

from <5 to IR. 

(3) For g G V* and i/ G Ad, the control influence operator B„ defined by 

= ^[P v («,0](0 = g f s [ks(u,0 ]{t)dv{s) (3.2) 

satisfies G L 2 ((0. T); V * ). 

Remark 1. The kernel k s and control influence operator B t , defined in this manner can be 
readily incorporated in a framework for proving model well-posedness and the existence of 
a measure which minimizes an appropriate parameter estimation functional. Moreover, as 
illustrated in Figure 7, this kernel yields a technique for characterizing the nested hysteresis 
curves which are common in magnetostrictive applications. 

4 Characterization of Hysteresis for Structural Models 
with Magnetostrictive Transducers 

To illustrate the use of the operator framework described in the last section to characterize 
hysteresis in magnetostrictive materials, we consider the modeling of magnetostrictive trans- 
ducers mounted to a thick cantilever beam as depicted in Figure 8. The transducers are 
considered to be mounted to the clamps at the fixed edge of the beam so that mass load- 
ing from the actuators themselves is minimized. A rigid bar is used connect the end of the 
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Figure 8. Cantilever beam with magnetostrictive actuators. 


Terfenol-D rod in the transducer to the beam. By driving the transducers out-of-phase, bend- 
ing moments are generated in a manner which can be used to attenuate beam vibrations. As 
described in [12], this experimental setup has been used in initial experiments to determine 
the potential of magnetostrictive transducers as structural actuators. Due to limitations in 
models and control laws, driving currents in the experiments were restricted to a range in 
which linearized results could be employed. Even in this restricted regime, the results of [ 12 ] 
demonstrate the utility of the magnetostrictive transducers for structural applications. 

It should be noted that a beam is considered here both due to its previous experimental use 
and the relative simplicity of the resulting model. The operator techniques are quite general, 
however, and can be applied in a similar manner to structures comprised of plates or shells as 
well as structures undergoing large deformations which leads to nonlinear models (e.g., von 
Karman models). 

For modeling purposes, we take the beam to have length L width b and thickness /?. The 
density. Young's modulus, Kelvin- Voigt damping coefficient and air damping coefficient for the 
beam are denoted by pb,Ek,CD h and 7, respectively. The cross-sectional area of the Terfenol 
rod is denoted by A mag while the Young's modulus and damping coefficient for the Terfenol 
rod are denoted by E H The length and width of the connecting bar are denoted by ( r 
and 6 r , respectively, while the bar density is given by p r . 


4.1 Strong Form of Beam Model 

The Timoshenko equations 






dM 


mag 


dx 


{t,x) 


w(t, 0) = a(t, 0) = 0 


t > 0 


0 < 0- < t , t > 0 

( 4 . 1 ) 


where to and a denote the transverse displacement and cross-sectional rotation, respectively, 
are used to model the beam dynamics. Exogenous surface forces to the beam are denoted by 



f{t,x) while p(x) is the composite density of the structure. The cross-sectional area of the 
beam is given by /4(;r ), while r(x) = I(x)/A(x) where I(x) is the moment of inertia of the cross- 
sectional area. Note that these quantities are spatially variable due to the nonhomogeneity 
in the region of the connection rod. Finally, the internal bending moment and shear force are 
given by M(t,x) arid Q(f,;r), respectively, while the external bending moment generated by 
the magnetostrictive actuators is denoted by M mag (t,x). 

To determine appropriate functional forms for the density, internal moment and shear force, 
the structural contributions due to the connecting bar and Terfenol rod must be quantified. 
We assume here that the connecting bars are perfectly rigid and contribute mass to the beam 
but do not affect the bending moments (we neglect air resistance to the bars). The actuator 
and Terfenol rod are considered to be supported from the boundary clamps so they do not 
contribute mass to the beam. The Terfenol rod is assumed to contribute an elastic stress 
which is uniform across the cross-sectional area of the rod. 

Under the assumption of uniform cross-sectional strains in the magnetostrictive rods, the 
density, stiffness and Kelvin- Voigt damping parameters for the structure are then taken to be 


p{ x ) — p^hb -(- 2p r b r t r \ ro d ( 3 ) 


E h 3 b 

EI(x) — — — b 2A mag E H (h/2 + t r ) 2 Xrod(x) (4*2) 

c D I(x) = C -~. + 2A mag Cp ( h/2 + t r f Xrod{x) 

where the location of the rods is delineated by the characteristic function \ ro d which has a 
value of 1 in the region covered by the connection bar and is 0 elsewhere. The internal moment 
and shear are then given by 


do d^o 

M(t,x) = EI(x)—{t,x) + c D I(x )— — (t,x) 


dx 


' dxdv 


c)3 

Q{t,x) = KAG(x)i3(t,X ) + kAcq{x) — {t,x) 


(4.3) 


where k is a correction factor which accounts for the fact that the outer surface of the beam 
cannot support a shear stress, G(x) is the shear modulus and cq(j) represents resistance to 
the shear strain rate. Finally, the shear deformations are defined by 


, dw , 

(3(t,x) = — — (t.x) + a(t,x) . 
dx 


It should be noted that the contributions due to the connection rods are dependent upon 
the exact experimental setup and different assumptions and models can also be used to in- 
corporate the passive rod contributions. In all cases, the piecewise constant parameters 
p(x), r(.r), EI{x), G(x), cpl(x), c.q(x) and constant parameter 7 must be estimated through 
a least squares fit to experimental data to attain a reasonable model for the specific experi- 
mental device (The values determined by (4.2) cannot be used with certainty when modeling 
the experimental apparatus due to inaccuracies in manufacturer specifications, etc.; however, 
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they can be used as initial values for the optimization routine.) Note that in the moment and 
shear expressions (4.3) and equations (4.1), parameters have been combined to yield a single 
value to be estimated. 

To characterize the external moment generated by the magnetostrictive transducers, we 
will assume that an offset DC current 7 0 is used to provide the magnetic field bias necessary 
to attain bidirectional strains. Ampere's law then yields the magnetic induction 

B{t) = fi*n[I(t) + Io] (4.4) 

where again, n is the number of terms per unit length in the solenoid, 1(f) is the current 
applied to the solenoid, and [i* is the permeability. As noted in Section 2, the permeability 
exhibits both stress-dependent nonlinearities and hysteresis due to the magnetic field. 

The strain can be related to the magnetic induction through the nonlinear equation 

e(t) = K B B(t) (4.5) 

where K B depends upon the magnetic induction B. Moreover, the external stress-strain 
relation 

<r(0 = E H e(t) (4.6) 

is also nonlinear since t he Young s modulus E n for the Terfenol-D rod can be highly dependent 
upon the applied magnetic field. Finally, the external moment generated bv the Terfenol-D 
rod is 

M mag (t,x) = 2 A maa (£ r + h/2)a(t)\ rod {x) (4.7) 

where f r denotes the length of the connection bar and A mag is the cross-sectional area of the 
Terfenol rod. 

Combination of (4.7) with (4.4)-(4.6) will yield a relationship for the external moment in 
terms of the applied current I(t). However, this relationship is not directly useful for modeling 
transducer dynamics since the parameters ft a ,k B arid E h are highly nonlinear and exhibit 
significant hysteresis. 

An alternative is to characterize the external moment through the fitting of hysteresis and 
material nonlinearities in terms of the Preisach kernels described in Section 3. Specifically, 
the external moment generated by the magnetostrictive transducer can be described by the 
relation 

Mmag(t,x) = 2 A mag (( r + h/2)[P „( /, £)](< )\w(*) (4.8) 

where 

[P„(J,O](0 = J s [k.(I,mt)Ms)- (4.9) 

The kernel k s is the extension of that defined in (3.1) to / £ CfO,! 1 ]. In essence, the method 
characterizes the hysteresis and material nonlinearities through curve fitting in terms of a 
measure v which, for a given set of experimental operating conditions, is estimated through 
a least squares fit to data. Due to its generality, the technique can be used to characterize 
a wide variety of transducer responses. Furthermore, the technique provides a framework 
suitable for parameter estimation and subsequent controller design. 
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4.2 Weak Form of Beam Model 


As noted when (4.3) and (4.8) are employed in (4.1), the use of the strong form of the beam 
model leads to the differentiation of discontinuous material parameters and inputs. It also 
necessitates the use of high-order approximating elements. To alleviate these difficulties and 
provide a framework amenable to analysis and approximation, we consider a corresponding 
weak form of the model. 

For this system, the state is taken to be y — (re, a) in the Hilbert space H — L 2 (0,/') x 
L 2 (0,() with the inner product 


= I p<j> x t x dx + [ prfafadx 
J o Jo 


where $ = {<i>\ , (p 2 ) , ® — (01 ,^ 2 )* The space of test functions is V — x Hj(Oj’), 

Hl( 0,f) = {4> £ H l (0,() | 0(0) = 0}, with the inner product 


($1 9)v — (<f>[ + 4 > 2)( 1 l ) l + V 7 2 ) dx + 0202 ■ 


It should be noted that with these choices, V is continuously and densely embedded in H and 
that V and H form a Gelfand triple; that is 


v v * . 


A weak form of the modeling equations is then 

f pwtfri dx + [ 7 wtl'i dx + f Qt[ dx + / prail '2 dx + / M 4 ' 2 dx + [ Qx [ } 2 dx 

Jo Jo Jo Jo Jo Jo 

= f M magi' 2 dx + / /01 dx 
Jo Jo 


(4.10) 


for all ^ = (yq, 4' 2 ) € V. In this form, derivatives are transferred onto suitably smooth 
test functions. This alleviates the difficulties associated with the discontinuities and reduces 
smoothness requirements on approximate solutions. 


4.3 Model Well-Posedness 

To provide a framework amenable to proving the well-posedness of the model and parameter 
estimation problem, it is advantageous to formulate the model in terms of sesquilinear forms 
and the operators which they generate. To this end, recall that y(t) = (w(t, •), a(t, *)) and 
define 


a x (q)(y{t)^) = f EIai^dx + 
Jo 

[q){y(t)iW)= f c D Ia f 44 2 dx + 
Jo 


f kAG{w ' + a)( 0j + 02 ) dx 
Jo 

I kAcq(w ' + a)(0i + i' 2 )dx + / 
Jo Jo 


dx 


(4.11) 
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where q = ( p , r, El , G, CpI, eg, 7) is considered in an admissible parameter space Q. It, can be 
directly verified that o\ and a 2 satisfy the boundedness and ellipticity conditions 

(HI) < Cl |$|v|«|v , c, € IR 

(H2) Recr^)^, $) > c 2 |$|v , c 2 > 0 


(H3) , C 3 € IR 

(H4) Recr 2 (<7)($,^) + c 4 \$\ 2 h > c 5 \$\y , c 4 6 1R,c 5 > 0 

for all <I>, G V where the constants rj, • • ■ , C 5 are independent of the parameters q G Q< 
From (4.8), it follows that the transducer contributions can be represented in terms of the 
operator B u : V — > V * by 

([B„(u,mt)i*)v.v = {[P»(v,t)](t)9,*)v,v 

= [P„(u,f)](0 f 2 A mag ( ( r + h j‘2 ) V’ 2 \ dx 

Jo 

where U — IR,, ('-, m )y*y is the usual duality product and P l is defined in (4.9). Note that 
from Theorem 1 , it follows that Z?„(t£,£) G L 2 ((0,T);V m ). Finally, with the definition /(f) = 
(p _1 /(/, *),0), we can write the weak form (4.10) in the abstract variational form 

V) v . v + ff a(g)(y(<)i *) + °i(q)(y(t), + f{t), V) v , v 

(4.12) 

y(0) = j/o , y(0) = yi 


for all T G V. 

An equivalent system can be obtained by invoking the boundedness of < 71,02 to define 
operators A t (q) G C{\\ V*) , i = 1,2, by 

<A,^») v . v = (T 1 (ijr)($,$) 

, 9 € V . 

(A 2 *,9) v . tV = <T 2 (q)(*,*) 

In operator form, the equations governing the beam dynamics then have the form 


VO) + Mq)y(t) + Mq)y(t) = [£„(«, f)](0 + /(<) (4.13) 


in V *. 

The subsequent result concerning the existence, uniqueness and regularity of solutions 
follows directly from Theorem 2.1 and Remark 2.1 of [4]. 


Theorem 2. Let Q denote a compact metric space and A4 denote the set of finite, signed 
Borel measures on 6*. Consider inputs u G C[0,T] and / G #(<S,{ — 1, 1}). Finally, let ai,a 2 
be given by (4.11) and hence satisfy (H1)-(H4). For each (</, v) G Q x A4, there then exists a 
unique solution y to (4.12) which satisfies 

y € C'((0, 7 7 ); V”) C /, 2 ((0, T); V) 

yeC(( 0.T);//)nL 2 ((0,7');V) 

yeL 2 ((0,T):V). 
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Remark 2. As noted in Remark 4.2 of [ 8 ], the regularity results in Theorem 2 can actually 
be strengthened to yield 

yeC(( o,nv) 

y € C((0, T); H) . 


4.4 Parameter Estimation 


The goal in the parameter estimation problem is to determine material parameters q = 
(p, r, El, G, cpl, cq, 7 ) and a measure v given data measurements z from some observable 
subspace Z of the state space. The form of the quadratic functional to be minimized depends 
on the experimental data which is available. For time domain data consisting of position, 
velocity or acceleration measurements at points x on the beam, an appropriate functional is 


J(q,v) 


E 


V s y 

dt s 


2 


(ti,x;q,v) - ii 


(4.14) 


with $ = 0, 1 or 2. Other functionals incorporating either time or frequency domain data can 
be found in [ 8 ]. In each case, the minimization is performed subject to y satisfying (4.12) for 
q £ Q and v £ A C A4 where A is an appropriate class of measures. 

As noted previously, the parameters p,r,EI ,G, cdI ,cq are assumed to be piecewise con- 
stant with partition points at the connection rod edges. The air constant 7 is constant, and 
all seven parameters are positive. The admissible parameter space Q is then taken to be a 
compact subset of the metric space Q = [A°°(0, Z ")] 6 x IR with elements piecewise constant 
between partition points and satisfying the positivity constraints. 

An appropriate choice for A is a set of probability measures with a metric which yields 
convergence in distribution. Specifically, let <Sa be a compact subset of S and let V{S a) 
denote the set of Borel probability measures v on <Sa* The space A is then defined to be 
V(S A ) endowed with the Prohorov metric which is defined for v \,u 2 € V(S±) by 


p{v\,V 2 ) = inf > 0 ( V\(F) < V 2 (F £ ) + £, F closed, F C S A ] 


(see [ 10 ] for details regarding this metric). Here F* denotes an s neighborhood for F. 


Remark 3. As detailed in [3, 6, 7, 10], the space A satisfies the following properties. 


(a) Let <7($a) denote the space of continuous functions on <Sa- Since V(S A ) C [C(«Sa)]*, 
convergence in the Prohorov metric is equivalent to weak* convergence; that is 


i'k v / fdu k -> / fdv 

Js A Js A 


for / G C{S A ). 


(b) A is a compact metric space. 
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(c) Consider the control influence operator B y defined in (3.2) and let y(q, v) be a solution to 
(4.12). For the functional J(q,v) defined in (4.14) witli fixed q £ Q, the map v J(q, v) 
is weak* lower semicontinuous from V(S&) to IR. 

The following theorem taken from [6] specifies conditions under which the parameter esti- 
mation problem with the magnetostrictive material inputs is well-posed. When combined with 
theory from [5], it provides a, framework for numerically estimating physical and hysteresis 
shape parameters through a least squares fit to experimental data. 

Theorem 3. Let B lf denote the magnetostrictive input operator defined in (3.2) and y(q, v) 
be the solution to (4.12). If <Sa is a compact subset of <5, then there is a probability measure 
v o G V{S a) which solves the minimization problem 

,/{q,u 0 )= inf J{q,v) 

(<Sa ) 

for the functional J{q,v) defined in (4.14). 


5 Concluding Remarks 

This paper addresses the characterization of material nonlinearities and hysteresis inherent to 
magnetostrictive materials at middle to high range drive levels. The hysteresis is induced by 
both magnetic fields and stresses within the material while additional nonlinearities arise in 
the strain-magnetic induction relations and stress-strain relations. Some of the mechanisms 
leading to these nonlinear input/output relations have been modeled in terms of electromag- 
netic and magnetomechanical theories [20, 21, 22]. However, several mechanisms governing 
magnetostrictive properties are still not. completely characterized by physics-based models. 
This motivates a mathematical characterization in terms of phenomenological Preisach tech- 
niques. 

Preisach models are empirical in the sense that they can be used to mathematically rep- 
resent hysteresis curves in terms of shape parameters determined through least squares fits 
to experimental data. Since the models are not directly derived from the physics of the sys- 
tem, they can be used to model dynamics in which the underlying physics is not thoroughly 
understood. The price paid for this generality is the loss of insight which can sometimes be 
provided by a physic-based model. 

In this work, Krasnoselskii-Pokrovskii kernels are used to characterize magnetostrictive 
inputs to a flexible structure. As detailed in [6], such kernels are advantageous over classical 
Preisach kernels due to the property that they are continuous with respect to both time 
and shape parameters. The operators generated via these kernels are then used to formulate 
the models in a Hilbert space framework. Well-posedness of the models is obtained through 
the theory of [4] while existence of optimal physical parameters and an optimal measure for 
characterizing the hysteresis curve follows from the theory of [3, 5, 6]. Hence for structural 
systems with inagnetostrictive transducers, this characterization technique provides a model 
which is amenable to analysis and approximation. In future investigations, the practical 
efficacy of this model will be tested through validation experiments. 
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